201903.00080v1 


chinaXiv 


人 Ab 


上 贱 


李文华 ” 男 1973 自 
生 ,， 教授 , ee 
器 可 靠 性 及 检测 技术 、 


ChinaXiv 合 作 期 刊 


航天 继电器 多 参数 贮存 寿命 加 速 退化 


水 田 


沈 


能 源 发 电 和 网 技术 


全 


培根 “ 男 1990 年 


， 硕 


研究 生 ， 研 究 方 


术 。 


向 为 电器 可 靠 性 及 检测 技 


预测 方法 研究 
李文华 沈 培 根 马 思 宁 
(河北 工业 大 学 电气 工程 学 院 天津 ”300130) 


摘要 : 航天 继电器 在 国防 、 航 天 工业 等 重要 领域 使 用 范围 十 分 广泛 ， 因 此 ， 为 了 
提高 航天 继电器 的 可 靠 性 ， 研 究 其 贮存 寿命 具有 十 分 重要 的 意义 。 目 前 ， 基 于 性 能 退 
化 数据 的 大 多 文献 局 限于 研究 单个 特征 参数 的 情况 ， 然 而 许多 产品 具有 多 个 性 能 退化 
特征 量 。 针 对 航天 继电器 具有 多 个 性 能 参数 的 特点 ， 本 文 提 出 了 退化 距离 分 析 法 ， 并 
在 此 基础 上 引入 主 元 分 析 法 (PCA)， 融合 了 多 个 特征 参数 ， 给 出 了 参数 权重 因子 求解 
方法 。 文 中 选择 了 航天 继电器 的 接触 电阻 RR、 吸 合 时 间 7 和 释放 电压 U, 这 三 个 参数 ， 
通过 该 方法 实现 了 贮存 寿命 预测 ， 并 选取 实测 数据 对 该 方法 进行 了 验证 。 
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Study on Multi-Parameters Storage Life Prediction Method 
for Aerospace Relay in Accelerated Degradation 


LiWenhua Shen Peigen Ma Sining 
(Hebei University of Technology Tianjin 300130 China) 


Abstract: Aerospace relay, is widely used in the field of national defense, aerospace 
industry, and others. Therefore, in order to improve the reliability of the aerospace 
relay, it is of great significance to study its storage life. At present, researches based on 
performance degradation data are mostly confined to the single characteristic parameters. 
However, many products have multiple performance degradation parameters. On account 
of the multiple performance parameters of aerospace relay, the degradation distance 
analysis method is proposed in this paper, principal component analysis (PCA) is also 
introduced into this analysis method to fuse multiple characteristic parameters, and 
weighting factors of parameters are given. In this paper, the aerospace relay contact 
resistance R;, pick-up time 7 release voltage Us are selected as life prediction parameters, 
and the relay storage life can be predicted with the help of the above method , which is 
also verified to be accurate by other data acquired in the experiment. 
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航天 继电器 广泛 用 于 国防 、 航 天 工业 等 重要 设 
备 中 ， 需 要 在 低温 、 高 温 、 高 温 及 高 频 振动 等 恶劣 
环境 下 完成 系统 电信 号 传递 、 电 路 隔离 及 继 电 保护 
等 重要 功能 。 在 实际 的 工作 情况 下 ， 只 要 有 一 台 继 
电器 出 现 故 障 ， 就 可 能 造成 整个 系统 瘫 病 ， 这 就 要 
求 继 电器 能 够 可 靠 工 作 以 保证 整个 系统 的 长 期 正常 
运行 。 有 些 继电器 的 贮存 时 间 远 远大 于 使 用 时 间 品 。 
在 贮存 过 程 中 ， 由 于 环境 条 件 的 不 断 变化 ， 继 电器 
的 性 能 会 出 现 衰退 ， 这 就 可 能 造成 当 其 投入 使 用 时 ， 
继电器 已 经 无 法 正常 工作 的 情况 。 因 此 ， 进 行 航天 
继电器 贮存 寿命 预测 对 保证 继电器 在 诸多 特殊 环境 
中 满足 极 高 的 可 靠 性 、 寿 命 要 求 和 消除 故障 隐患 可 
能 带 来 的 重大 的 经 济 损失 有 着 深远 的 意义 中。 传统 
的 寿命 试验 方法 对 继电器 进行 可 靠 性 评估 和 寿命 预 
测 时 ， 会 面临 试验 成 本 高 、 失 效 数 据 少 等 许多 问题 ， 
从 而 无 法 得 到 有 效 的 结果 。 继 电器 产品 在 服役 过 程 
中 的 失效 机 理 最 终 能 够 追溯 到 产品 潜在 的 性 能 退化 
过 程 ， 可 以 认为 性 能 退化 最 终 导 致 了 产品 失效 (或 
故障 ) 的 产生 。 因 此 ， 基 于 性 能 退化 数据 的 分 析 已 
经 成 为 一 种 有 效 且 节约 成 本 的 性 能 评估 方法 中。 
2 问题 描述 

产品 性 能 退化 可 能 会 引起 某 一 个 或 几 个 特征 性 
能 参数 发 生变 化 ， 这 就 需要 根据 产品 功能 来 考察 其 
中 某 几 个 重点 特征 参数 ， 然 后 综合 考虑 这 些 参数 来 
评价 产品 工作 状态 或 性 能 情况 。 产 品 的 性 能 特别 
是 一 些 复杂 的 系统 ， 往 往 要 几 个 性 能 参数 综合 反 
映 ， 同 样 需要 综合 考虑 这 些 参数 来 评价 产品 的 性 能 
退化 外 。 在 产品 性 能 退化 数据 分 析 研究 中 ， 针 对 多 
性 能 退化 参数 的 文献 大 多 假设 多 个 退化 参数 的 联合 
概率 密度 函数 服从 某 种 多 元 分 布 ， 例 如 文献 [3] 假 
设 检 查 时 刻 点 的 多 元 性 能 退化 数据 服从 多 元 正 态 
分 布 ， 文献 [5] 利用 了 联合 概率 密度 函数 表征 产品 
失效 概率 密度 ， 或 者 假设 性 能 参数 退化 规律 符合 
Wiener、 泊 松 等 随机 过 程 并 通过 预先 确定 的 单 特征 
参数 的 失效 阔 值 来 计算 产品 的 故障 时 间 “1 ， 这 就 忽 
略 了 多 性 能 退化 参数 相关 条 件 下 ， 失 效 阔 值 也 可 能 
变化 的 情况 。 本 文 提出 了 退化 距离 分 析 法 ， 解 决 了 
多 参数 和 多 参数 失效 阔 值 的 问题 。 在 利用 该 方法 进 
行 分 析 需 要 做 如 下 假设 : 

(1) 产品 的 各 个 性 能 参数 随时 间 退 化 具有 不 可 
逆 性 中 。 
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(2) 多 个 性 能 退化 参数 相关 或 者 相互 独立 。 
(3) 当 继电器 的 某 个 退化 参数 达到 产品 给 定 失 
效 国 值 时 ， 就 判定 产品 失效 。 


3 ”数据 预 处 理 


对 25 台 某 型 号 航天 继电器 在 125C 恒定 应 力 下 
进行 加 速 喧 存 退 化 试验 ， 采 集 继电器 的 接触 电阻 Rj、 
吸 合 时 间 未 和 释放 电压 U, 这 3 个 电 特 征 参数 ， 其 
中 选取 触 点 接触 电阻 R 的 最 大 值 为 产品 的 接触 电阻 
值 。 产 品 的 规定 国 值 分 别 为 30 mQ、8 ms、2.5 V。 
试验 中 编号 17# 继电器 贮存 试验 测量 得 到 的 特征 量 
数据 见 表 1。 


表 1 特征 参数 样本 数据 


Tab.1 Characteristic parameters of sample data 


周 数 接触 电阻 RimQ ” 吸 合 时 间 TWms ”释放 电压 UW/V 


1 8.31 5.09 4.50 
2 8.34 5.68 4.21 
8.66 $:82 4.13 
4 8.65 5.74 4.29 
3 8.23 $.76 4.15 
6 8.87 6.33 4.23 
7 8.56 5.81 4.10 
8 8&33 5.86 4.11 
9 8.62 3:82 4.09 
10 8.63 6.78 4.04 
11 8.93 $9 4.06 
12 8.66 3.92 4.13 
13 9.04 5.70 4.06 
14 8.70 $:85 4.06 
15 9.07 5.90 4.25 
16 8.86 SA 4.18 
ly 8.88 5.86 4.33 
18 8.91 5.96 4.13 
19 8.74 5.85 4.28 
20 9.08 6.02 4.21 
21 9.00 3.9] 4.23 
22 8.99 5.91 4.05 
23 9.20 5.82 3.96 
24 9.50 95 3:97 


由 于 测量 设备 本 身 存 在 误差 和 外 界 环境 等 扰动 
因素 ， 所 测 的 特征 变量 包含 一 定 的 噪声 ， 为 减少 干 
扰 的 影响 ， 降 低 测量 值 的 随机 性 ， 前 期 有 必要 对 数 
据 进行 平稳 处 理 ， 本 文采 用 移动 平滑 滤波 法 分 别 对 3 
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个 特征 参数 进行 平滑 处 理 。 
4 ”参数 退化 距离 分 析 


在 大 多 数 实际 问题 中 ， 表 示 产 品 工作 状态 的 性 
能 参数 ， 其 稳定 的 变化 趋势 是 单调 上 升 或 单调 下 降 
的 。 从 产品 性 能 开始 发 生 退 化 到 产品 最 终 失 效 ， 产 
品 的 性 能 参数 初 值 是 朝 着 失效 国 值 方向 发 展 的 。 随 
着 试验 时 间 的 延长 ， 每 个 性 能 参数 值 与 失效 国 值 间 
的 距离 将 越 来 越 接 近 ， 当 参数 值 无 限 接近 失效 国 值 
的 时 候 ， 也 就 是 距离 接近 为 0 时 ， 产 品 最 终 发 生 了 
失效 。 这 种 现象 也 反映 了 退化 过 程 的 不 可 逆 性 。 用 
2 个 性 能 参数 全 ，X 来 说 明 产 品 性 能 退化 过 程 ， 如 
图 1 所 示 。 


性 能 参数 初始 值 ( 双 ， 矶 ) 
图 1 产品 性 能 退化 过 程 


Fig.l Product performance degradation process 


假设 产品 有 ZL 个 性 能 退化 参数 ; X', XX， 
XY ，…X“ 且 设 第 j 个 性 能 退化 参数 的 失效 六 值 为 
守 ， 则 在 上 时 刻 性 能 退化 参数 值 与 失效 闽 值 间距 离 
为 芝 (0- 嫩 ， 记 做 DOD， 在 上 时刻 多 性 能 参数 与 效 
值 间 距离 的 定义 为 


p=- Sw (1) 


式 中 ，w/(j=1, 2,…Z) 为 加 权 因 子 ， 根 据 各 个 特征 
参数 在 实际 中 的 重要 程度 给 出 。 通 过 对 某 时 刻 的 参 
数值 与 失效 国 值 间距 离 的 分 析 ， 可 以 看 出 距离 可 以 
有 效 地 反映 产品 性 能 的 变化 趋势 ， 那 么 在 实际 问题 
中 ， 如 何 从 若干 特征 参数 中 筛选 较为 重要 的 参数 ， 
确定 参数 之 间 的 加 权 因 子 和 经 数据 变化 后 失效 阔 值 ， 
是 进行 性 能 退化 分 析 的 关键 所 在 。 
4.1 主 元 分 析 融 合 特征 参数 

主 元 分 析 是 一 种 常用 的 特征 提取 方法 ， 已 广泛 
应 用 于 传感器 信号 的 提取 、 发 电量 预测 和 凝 汽 器 故 
障 诊断 等 领域 。 采 用 主 元 分 析 法 能 把 多 维 的 特征 参 
数 维度 降低 ， 能 提取 出 所 有 特征 参数 之 间 的 相关 性 
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较 强 的 特征 量 。 其 基本 思想 是 通过 对 原始 数据 进行 
线性 变换 ， 变 换 后 的 数据 即 主 元 ， 上 共有 较 低 的 维 数 
但 能 最 多 地 反映 原始 数据 信息 " "“。 对 于 某 个 数据 
和 矩阵 闷 xw， 每 一 列 天 对 应 一 个 特征 矢量 ， 每 一 行为 
观测 样本 。 对 和 矩阵 也, 进行 主 元 分 析 实 质 上 就 是 对 
其 协 方差 矩阵 进行 特征 矢量 分 析 。 当 马 *w 的 特征 矢 
量 量 纲 不 同 或 数量 级 差别 较 大 时 ， 需 采用 式 (2) 对 
其 进行 数据 极 差 变换 ， 经 过 变换 后 ， 和 矩阵 关中 的 每 
个 元 素 取 值 均 在 0 ~ 1 之 间 。 


Xi, 一 InaXX 

起， 7 1<k<n 
JU = 一 一 一 一 
max Xx,; — min Xi 
I<k<n YY 1<k<n 了 


订 


(2) 


本 文 以 接触 电阻 R、 吸 合 时 间 和 释放 电压 
及， 这 3 个 特征 参数 构成 数据 矩阵 Xiyx3。 对 Xx， 
进行 主 元 分 析 ， 得 到 各 主 元 Y(i= 1 2, 3)。 当 各 特征 
参数 之 间 存 在 一 定 的 相关 性 时 ，X; 的 变化 主要 体 
现在 前 儿 个 主 元 上 。 经 PCA 计算 得 到 3 个 主 元 的 系 
数 和 矩阵 


0.0444 0.0735 0.9992 
0.8702 0.4868 -0.0757 
-0.4907 0.8702 -0.0441 


在 应 用 主 元 分 析 法 进行 数据 压缩 和 特征 提取 时 ， 
为 了 使 舍弃 的 原 有 数据 信息 不 影响 对 整体 数据 的 分 
析 ， 要 求 主 元 的 累积 页 献 率 必须 大 于 85%， 见 表 2。 

表 2 主 元 贡献 率 


Tab.2 Principal component contribution 


主 元 特征 根 贡献 率 (%) 
马 1.10x 107 97.52 
了 2.57 x 10™ 228 
2.16 x 107 0.20 


第 一 主 元 和 第 二 主 元 的 累积 贡献 率 已 经 达到 
85% 以 上 。 经 过 PCA， 不 但 消除 了 多 参数 间 彼 此 的 
元 余 信 息 及 干扰 信息 ， 而 且 第 一 主 元 和 第 二 主 元 已 
经 足够 强 地 反映 出 原始 多 维特 征 参数 的 变化 信息 。 
所 以 ， 完 全 可 以 省 略 第 三 主 元 ， 只 保留 第 一 、 第 二 
主 元 ， 从 而 达到 降低 多 参数 维度 的 效果 。 第 一 主 元 
的 贡献 率 最 大 ,保持 了 特性 参数 之 间 的 绝 大 部 分 相 
关 性 ， 反 映 了 特征 参数 的 变化 。 在 忽略 其 他 贡献 率 
较 小 主 元 的 情况 下 ， 第 二 主 元 是 原始 信息 和 第 一 主 
元 之 间 误 差 信 息 的 一 种 综合 体现 ， 反 映 了 参数 变化 
中 不 确定 性 和 干扰 因素 。 由 系数 矩阵 ， 可 得 到 第 一 
主 元 五、 第 二 主 元 也 的 表达 式 
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¥, = 0.0444x +0.8702x, -0.4907x; (3) 


,=0.0755x, +0.4868x, +0.8702x, (4) 


根据 式 (3) 和 式 (4) 可 以 进一步 得 到 第 一 主 
元 和 第 二 主 元 变化 的 时 间 序 列 。 
4.2 主 元 参数 权 因子 的 确定 方法 

第 / 个 主 元 的 贡献 率 代表 的 含义 是 第 /个 主 元 
所 包含 信息 占 系统 全 部 信息 的 百分比 ， 即 第 /个 主 
元 的 重要 程度 的 百分比 。w 表示 式 (1) 中 ， 各 个 
参数 的 重要 程度 ， 也 就 是 占有 信息 量 的 多 少 。 若 先 
取 了 7 个 主 元 ,可 得 


oo- Sn (5) 


式 中 , 4 表 示 为 被 选取 的 第 j 个 主 元 的 特征 值 ， 这 
里 提取 了 两 个 主 元 了、 钱 。 根 据 式 (5)， 可 以 计算 
出 参数 权 因 子 w,、w; 的 数值 。 
4.3 多 参数 失效 阅 值 的 确定 

首先 ， 在 分 析 性 能 退化 的 数据 时 ， 为 了 消除 
各 个 性 能 参数 之 间 量 纲 的 不 同 对 原始 数据 进行 了 
变换 。 其 次 ， 对 消除 量 纲 后 的 数据 进行 了 主 成 分 
分 析 。 因 此 ， 在 给 定 各 个 性 能 参数 失效 阅 值 的 情 
况 下 ， 不 能 简单 地 得 到 失效 国 值 ， 而 是 需要 对 失效 
国 值 也 做 上 述 相应 的 变化 ， 才 能 保证 失效 浆 值 的 有 
效 性 。 

在 产品 性 能 发 生 退 化 的 过 程 中 ， 从 长 期 来 看 ， 
产品 的 性 能 指标 值 的 变化 趋势 是 上 升 或 下 降 的 ， 而 
这 种 变化 趋势 就 是 产品 的 性 能 退化 轨迹 ， 显 然 产品 
的 各 个 参数 的 失效 国 值 就 是 退化 轨迹 上 的 最 大 值 或 
最 小 值 。 根 据 极 差 变换 公式 和 数据 所 反映 的 趋势 ， 
可 以 对 产品 说 明 给 定 的 失效 靖 值 进行 变换 。 接 触电 
阻 Rj、 吸 合 时 间 7 和 释放 电压 VU,，3 个 特征 参数 规 
定 阔 值 向 量 为 [50 mQ, 8 ms, 2.5 V]， 则 经 过 极 差 变 
换 后 国 值 向 量 为 [1, 1, 0]。 通 过 PCA 分 析 ， 利 用 式 
(3) 和 式 (4) 把 3 个 特征 参数 的 初始 国 值 变 为 两 个 
主 元 的 国 值 。 得 出 第 一 主 元 国 值 Y= 0.91， 第 二 主 
元 国 值 及,= 0.56。 

4.4 主 元 退化 距离 分 析 

将 得 到 第 一 、 二 主 元 时 间 序 列 蕊 (D、 到 (0D) 各 点 
值 和 由 式 (5) 得 到 的 权重 因子 带 入 式 (1) 中 得 到 
的 数值 见 表 3。 

对 表 中 各 点 进行 描 点 ， 并 用 最 小 二 乘法 进行 曲 
线 拟 合 ， 如 图 2 所 示 。 


ChinaXiv 合 作 期 刊 


表 3 主 元 参数 退化 距离 的 时 间 序 列 
Tab.3 The time series of principal parameter 


degradation distance 


周 数 主 元 参数 距离 周 数 主 元 参数 距离 
1 1.373 13 1.052 
2 1.192 14 1.073 
3 1.161 15 1.087 
4 1.076 16 1.075 
5 1.063 17 1.086 
6 1.059 18 1.077 
7 1.045 19 1.068 
8 0.981 20 1.051 
9 0.993 21 1.051 
10 0.988 22 0.949 
11 0.994 23 0.881 
12 0.991 24 0.850 

+ 散 点 
一 拟 合 


主 元 参数 距离 
Se 一 二 一 一 一 一 一 一 
己 一 已 上 个 ~10o 习 品 一 忆 wm 上 内 
人 


人 
un 


10 13 20 25 30 
周 数 

2 最 小 二 乘法 拟 合 
Fig.2 Least square fitting 


可 以 看 出 离散 点 均匀 的 分 布 在 一 条 直线 两 侧 ， 
根据 曲线 的 趋势 利用 线性 方程 对 各 个 离散 点 进行 拟 
合 ， 得 到 的 拟 合 函数 为 


D(1)= pit+p; (6) 


式 中 ， 函 数 系数 为 : 站 =-0.008 867 9, p,=1.161 5。 

当 产品 发 生 失 效 时 ， 主 元 参数 国 值 距离 函数 
DOD=0。 可 以 算出 4=130， 即 在 125C 的 应 力 条 件 
下 该 产品 伪 失 效 寿 命 为 130 周 。 在 实际 应 用 中 ， 相 
同 试验 条 件 下 一 批 同型 号 的 继电器 寿命 差距 是 不 大 
的 ， 所 以 本 文 的 思路 是 把 该 方法 用 到 该 试验 条 件 下 
其 他 航天 继电器 退化 数据 上 进行 类 比 验 证 ， 如 果 在 
其 他 继电器 上 能 得 到 相差 不 大 的 寿命 值 ， 则 说 明 该 
方法 具有 一 定 适 用 性 和 准确 性 。 从 试验 的 25 台 继 
电器 中 随机 选择 编号 为 3#、5#、8#、13# 这 4 台 继 
电器 ， 以 所 测 得 数据 为 例 对 该 方法 进行 类 比 验证 ， 
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得 到 这 4 台 继 电器 主 元 参数 退化 距离 的 时 间 序 列 
见 表 4 D 
表 4 4 台 主 元 参数 退化 距离 的 时 间 序 列 
Tab.4 The time series of principal parameter 


degradation distance 


主 元 参数 距离 
周 数 

3# 5# 8# 13# 
1 1.367 1.451 1.453 1.324 
2 L212 1.340 1.366 1.240 
3 1.135 1.224 1.336 1.212 
4 1.097 L137 1.288 1.166 
5 1.094 1.134 1.272 1.151 
6 1.092 1.134 1.259 1.141 
7 1.091 :139 1253 1.137 
9 1.075 1:129 1.221 1.109 
0 1.067 1.117 1.214 1.105 
11 1.054 1.105 1.216 1.108 
2 1.046 1.098 1.203 1.098 
3 1.045 1.094 1.221 1 
4 1.030 1.092 1.220 1.119 
入 1.024 1.088 1.214 1 .13 
16 1.010 1.082 1.198 1.101 
了 1.008 1.074 1.192 1.098 
8 1.006 1.067 1.186 1.094 
9 0.996 1.063 1.176 1.086 
20 0.994 1.063 1:170 1.082 
21 0.991 1.062 1.163 1.077 
22 0.982 1.052 1.149 1.065 
23 0.978 1.030 L136 1.054 
24 0.988 1.018 1.114 1.033 


根据 表 4 所 示 的 4 台 主 元 参数 退化 距离 时 间 序 
列 进行 了 寿命 预测 ， 得 到 3#、5#、8#、13# 继 电器 
伪 寿 命 值 分 别 t=117、 t=114、 t=140、 t=150。 
设 由 表 3 所 求 17# 继 电器 的 伪 寿 命 值 为 t， 则 5 台 
继电器 伪 寿 命 均值 T=(t+b+h+ty+1s)/5 =130.2，; 
伪 寿 命 标准 差 : D=std(t)= 15.2。 由 均值 和 方差 可 以 
看 出 ， 将 125 的 应 力 条 件 下 随机 抽取 的 5 台 继 电 
器 数据 带 入 该 方法 所 预测 出 的 寿命 值 相差 不 大 。 这 
表明 该 方法 用 到 其 他 航天 继电器 退化 数据 上 进行 类 
比 验证 时 ， 可 以 得 到 相差 不 大 的 寿命 值 ， 说 明 该 方 
法 能 够 适用 于 其 他 继电器 数据 ， 并 具有 一 定 准 确 性 ， 
能 够 实现 基于 多 参数 航天 继电器 性 能 退化 的 贮存 寿 
命 预 测 。 
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5 结论 


(1) 加 速 贮 存 试验 能 够 在 相对 短 时 间 内 获得 关 
于 继电器 工作 特性 退化 的 关键 信息 ， 通 过 对 退化 数 
据 的 分 析 ， 能 够 完成 航空 航天 继电器 的 贮存 寿命 
预测 。 

(2) 本 文 提供 了 一 种 分 析 多 性 能 参数 退化 数据 
的 新 思路 ， 即 主 元 参数 退化 距离 分 析 法 。 给 出 了 参 
数 与 相应 国 值 间 退化 距离 的 定义 ， 并 引入 了 主 元 分 
析 法 ， 给 出 了 不 同 主 元 参数 在 退化 距离 分 析 中 加 权 
因子 的 确定 方法 。 最 后 ， 预 测 了 航天 继电器 的 贮存 
寿命 。 

(3) 解决 了 现 有 方法 基于 多 元 联合 概率 密度 函 
数 假设 ， 以 及 求解 过 程 较 复杂 等 不 足 。 特 别 是 当 产 
品 的 性 能 参数 增多 时 ， 主 元 参数 退化 距离 分 析 法 能 
大 大 减 小 计算 量 。 

(4) 虽然 线性 拟 合 的 方法 能 够 直接 快速 地 通过 
主 元 退化 距离 预测 出 产品 的 妨 存 寿命 ， 但 是 拟 合 的 
好 坏 直接 依赖 距离 数据 的 质量 和 数量 ， 存 在 一 定 的 
局 限 性 。 因 此 ， 还 有 待 进一步 研究 和 完善 。 
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